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Аннотация 

Введение. Рассматривается задача о критических нагрузках сжатой ортотропной прямоугольной пластины на 
упругом основании. Для пластины заданы параметры ортотропии: коэффициенты Пуассона, модули Юнга для 
главных направлений и модуль сдвига материала пластины. Составляющие сжимающей нагрузки равномерно 
распределены по двум противоположным краям пластины и действуют параллельно осям координат. Края 
пластины свободно защемлены или шарнирно оперты. Рассмотрены также случаи, когда два параллельных края 
пластины свободны от нагрузок, а два других свободно защемлены или шарнирно оперты. 

Материалы и методы. Задача рассматривается на основе системы нелинейных уравнений равновесия типа 
Кармана. Критические значения параметра нагрузки определяются из линеаризованной на тривиальном 
решении задачи. При этом для решения краевой задачи на собственные значения применяется вариационный 
метод в сочетании с конечно-разностным методом. 

Результаты исследования. Задача сведена к решению параметрической линейной краевой задачи на 
собственные значения. В случае краевых условий подвижного шарнирного опирания приведены точные 
формулы собственных значений и собственных функций, а в случае свободного защемления краев применен 
вариационный метод в сочетании с конечно-разностным методом и построена компьютерная программа 
решения задачи. Установлено, что критическому значению параметра сжимающей нагрузки, при котором 
происходит потеря устойчивости сжатой пластины, может соответствовать одна или две собственные функции, 
выражающие прогиб пластины. Приведены результаты численных расчетов критических значений сжимающей 
нагрузки при различных значениях параметров ортотропии и построены графики соответствующих форм 
равновесия. Для случая длинной ортотропной пластины на упругом основании установлено, что главный член 
асимптотического разложения решения линейной задачи на собственные значения определяется из задачи о 
критических нагрузках сжатой балки на упругом основании с модулем упругости, совпадающим с модулем 
упругости пластины по продольному направлению. 

Обсуждение и заключения. Исследована задача о критических нагрузках сжатой в двух направлениях 
ортотропной пластины, лежащей на упругом основании. При повышении составляющей сжимающей нагрузки 
вдоль одного направления снижается величина критического значения нагрузки, сжимающей пластину вдоль 
другого направления. Если ортотропная пластина сжимается нагрузкой вдоль направления, которое 
соответствует большей изгибной жесткости, то критическое значение потери устойчивости больше, чем 
критическое значение действующей вдоль направления меньшей изгибной жесткости сжимающей нагрузки. 
Наличие упругого основания повышает несущую способность сжатой пластины. 


Ключевые слова: критическая нагрузка, упругая ортотропная пластина, устойчивость равновесия, упругое 
основание, параметрическая краевая задача на собственные значения. 
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Введение. Л. М. Зубов в работе [1] вывел уравнения равновесия упругой пластинки, содержащей источники 
внутренних напряжений в виде дислокаций и дисклинаций, являющиеся модификацией уравнений Кармана. В 
этой работе также решена задача об изгибе тонкой пластинки (мембраны) под действием внутренних 
напряжений, обусловленных дефектами. В работах [2—5] теория дислокаций и дисклинаций применяется в 
исследовании поведения пластин и оболочек. В статье [6] построено общее решение уравнения колебаний 
прямоугольной ортотропной пластины со свободными краями. В работах [7, 8] строятся асимптотические 
решения задач статики и динамики узких пластин. 
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В работе [9] рассмотрена задача о ветвлении равновесий под влиянием малого нормального давления 
сжатой в одном направлении прямоугольной пластины, содержащей дислокации и дисклинации. Рассмотрен 
случай, когда продольные края пластины свободны от нагрузок, а два других края закреплены. В [10] для 
задачи о критических нагрузках сжатой узкой пластины строится асимптотическое решение с помощью 
разложения по малому параметру относительной ширины пластины. В [11] исследована задача о ветвлении 
равновесий упругой прямоугольной пластины с внутренними источниками напряжений для случая, когда 
сжимающие нагрузки неравномерно распределены по краям. Установлено, что при четных формах функции 
несовместности и четных формах распределения по краям сжимающих усилий наличие малой нормальной 
нагрузки не снижает несущую способность пластины. 

В работе Морозова Н. Ф., Беляева А. К., Товстика П. Е., Товстик Т. П. [12] представлен асимптотический 
вывод двумерных уравнений равновесия тонкой упругой неоднородной пластины, изготовленной из 
анизотропного материала общего вида. В [13] сформулированы выводы о формах потери устойчивости сжатой 
пластины на упругом мягком основании. В работе [14] исследована жесткость на изгиб тонкой упругой 
многослойной пластины с трансверсально изотропными слоями. Проводится сравнение двух моделей учета 
влияния поперечного сдвига. 

В работе [15] построена модель термоупруговязкопластической деформации композита, который 
перекрестно армирован непрерывными волокнами в произвольных направлениях. 

В[6] в рамках прямого подхода к теории пластин рассматривается задача о гиперупругой пластине с 
неоднородно распределенными начальными напряжениями. Пластина рассматривается как материальная 
поверхность с пятью степенями свободы (три перемещения и два вращения). В работе [17] в рамках 
нелинейной упругости анализируется устойчивость равномерно сжатой круглой двухслойной пластины с 
первоначально сжатым или растянутым слоем. Для определяющего соотношения материала используется 
модель несжимаемого неогуковского материала. Представлен анализ зависимости результирующих 
критических напряжений от начальных деформаций и параметров жесткости. В работе [18] исследуется изгиб 
трехслойной пластины, содержащей поверхностные и межфазные напряжения, рассмотрена теория пластин с 
деформациями сдвига первого порядка и модель поверхностных напряжений Гуртина-Мердока. Анализируется 
зависимость параметров жесткости пластины от модулей поверхностной упругости. В работе [19] 
рассматривается задача о неустойчивости трехслойной нелинейно-упругой прямоугольной пластины с 
предварительно напряженным средним слоем. 

В настоящей работе исследуется влияние параметров материала пластины и параметра упругого основания 
на критические нагрузки потери устойчивости сжатой упругой прямоугольной ортотропной пластины на 
упругом основании. 

Материалы и методы. Рассмотрим лежащую на линейно упругом основании ортотропную прямоугольную 
пластину, которая сжимается усилиями Ри О вдоль осей Х и У. Уравнения равновесия можно записать в 
виде [7]: 

2,0%. \+20,0.0.\+р,04\М+КМ = [№, Е|, 


т 1 2м 
Е, сЕ 


(р 


1 1 
90 Е+ Е 0Е=- - №, М]. 


Краевые условия запишем в виде: 


У\У=д,\=0, 0,0,Е=0, дуЕ=-Р при |Х| =а/2, т 

У\У=д,\=0, 0,0.Е=0, дхЕ=-О при |У| =Ъ/2. 

У\У=д,\=0, 0,0,Е=0, д,Е=-Р при |[Х| =а/2, 8 

У\=д,\=0, 0,0,Е=0, дхЕ=-О при |У| =Ъ/2. 

\=д,\=0, 0,0,Е=0, дуЕ=-Р при |Х| =а/2, и 
УМ +0 М = 0, +(2-у,)0.0,\ =9,0,Е =дУЕ = 0 при [У| = 5/2. 

\=0.\=0, 0,0,Е=0, д,Е=-Р при |Х| =а/2, 8 
УМ +0, М = 9, +(2-у,,)0%0,\ =0,0,Е =0УЕ=0 при [У| =Ъ/2. 
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[\,Е] = 9х \0уЕ+0,\0хЕ-20,0,\0,0,Е, Г, = Е.В? /(12(1-ч,,)), 1=1,2; О, =©68*/2, О, =Буу, +20,, 
Ех, =Е м. 

Здесь у, \,, Е,, Е, С — коэффициенты Пуассона, модули Юнга для главных направлений и модуль 
сдвига материала пластины; О, О, — жесткости изгиба по главным направлениям; О, — жесткость 
кручения; функция У (Х, У) выражает прогиб, а Е (Х, У) обозначает функцию напряжений; й — толщина 
пластины; К\У/ — реакция основания; К — коэффициент постели; начало координат Х, У располагается в 
центре пластины и оси параллельны ее краям; 0х — частная производная и-го порядка по переменной Х. Мы 


полагаем, что края пластины свободно защемлены (2) или шарнирно оперты (3). При этом сжимающие усилия 
Ри О равномерно распределены по краям Х=+а/2 и У=+6/2 соответственно. Рассматриваются также 


случаи, когда два параллельных края пластины У =-6/2 свободны от нагрузок, а два других свободно 


защемлены (4) или шарнирно оперты (5). В этих случаях сжимающая нагрузка приложена только к краям 
Х=+а/2. 


2 2 
Положим Ва“ -0^ о Тогда для функции Ф справедливы однородные краевые условия 
2 2. 


Фу =Ф,, =0 при |Х|=а/2, Ф,=Ф,‚„ =0 при |У|=5/2. Эти условия можно заменить на эквивалентные 


краевые условия Ф =Ф, =0 при |Х|=а/2, Ф=Ф, =0 при |У =6/2. Перейдем к безразмерным переменным по 


формулам: 
Х=ах; У=Ьу, 6 ь ; Ф(Х,У) = Р.Е (х,у); \(Х,У) =м(х,у)в; = ; Р-Р ; 9-95 , 
а 
Е О, . р, 
а=6(1-му,); 2, = “2%: с, Е. м 
Тогда систему уравнений (1) и краевые условия (2), (3) можно переписать в виде: 
50 + 2с,6°0.0у\ + с,дум +5*Ку = 5* [м,Е]-р5?дхм - адм, © 
50+ 2.0.0.6" Е-+с,дуЕ = -ас,6* [м, м |. 
[№ де, Е, 9.1] „=0; [дум =! (7) 
‚дм, Е, 0 =0; | м, 0? \,Е, 0 =0. 
[№ "м, Е, 0, 1 [м (БО, . (8) 
а 2 2 Е 2 — 
[м,д, ч.1, 0,1], =0: [дм + м.б, дз + (2-м, ) 92 д,м,1, и], = (9) 
2 И 2 2 З И 2 ыы 
[№ 0? ч.1, 01 | диз 20; [дум дум +(2—м.)09:дум, 60] „= (10) 


При любых значениях параметров (, р, 4 нелинейной краевой задаче (6) с любым из условий (7)—(10) 
удовлетворяет тривиальное решение (м, Е.) = (0,0). 
2 .. ь и 
Обозначим через Е` множество вектор-функций Ё Е) с конечной нормой, которая определяется с 


помощью скалярного произведения: 


ее = [[(е2, +Ё,2, )ахду; Е аи (ыы (11) 


о 


Через Е' обозначим пространство бесконечно дифференцируемых в области @={(х, у): |‹|</2|У|< 1/2} 


вектор-функций и = (\,Е) ; у= (м,Е) ‚ которые удовлетворяют на границе ОО области © одному из краевых 
условий (5) или (6). Норму в Е! определим с помощью скалярного произведения: 
0 ду 
(альт, - (12) 
эк \0х'ду’ дх'ду' | 
Краевая задача (6) с одним из краевых условий (7)—(10) может быть рассмотрена как нелинейное 
операторное уравнение: 
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Ми =Па, и=(\м,Р)ЕЕ', (13) 
где 
5'0%\+2с,6°0.дум +с,дум + Км + бд +9дум А 5° [№] 


= - (14) 
5'0 +2. 00254,04 


Ю 


ой г 
—ас,5* [м,\] 


Здесь М, — линейный, а П — нелинейный операторы, так как: 
[№, Е] = деждЕ +д.м0Е—20,0, 0,04. 
Из результатов работ И. И. Воровича и Н. Ф. Морозова следует [9], что операторы М, и П действуют из 


пространства Е' в Е”. Пусть А = 0, р=0, 4 = 0. Тогда уравнение (9) имеет единственное тривиальное решение 


и, = (.,Е.) = (0,0). При малых изменениях величин параметров р, а и К тривиальное решение остается 
единственным. Если при любых малых значениях числа &>0 и фиксированном значении К существуют 
значения параметров р, 4, которые удовлетворяют неравенству |р-Р.|+|9-9|<= и для которых 
существуют хотя бы два решения уравнения (9), то точку (Ро, Чо; называют точкой ветвления (или точкой 


бифуркации) [9] этого уравнения. В работе [9] показано, что для нахождения точки бифуркации [у 4) надо 


решить краевую задачу на собственные значения, которая представляет собой линеаризованное на тривиальном 
решении уравнение: 
Ми =0. (15) 


Собственные значения (ро ›‚ Чо ) ЯВЛЯЮТСЯ критическими значениями параметров сжимающей 


нагрузки (р,а). Учитывая (13) и (14), уравнение (15) можно записать как следующую краевую задачу на 


собственные значения: найти такие значения параметра р, что для фиксированных значений параметра 
основания (коэффициента постели) Ки параметра 4 нагрузки вдоль оси у, следующая задача имеет 
нетривиальное решение: 


5‘ 0%\ + 2с,6°0.дум + с, дум + Км +адум = -рб*дхм, (16) 

[№.0, 1, =0; [6 ] =0, (17) 

[м,д,м в =0; [95%+ мдм, дум + (2,0: Е,01 |, =0, (19) 
[д | низ 20; [9% + ду, дум + (2, 0.0 ч.0] = (20) 


При этом критическому значению потери устойчивости сжатой пластины отвечает минимальное 
собственное значение параметра р, а соответствующая этому значению собственная функция показывает форму 
равновесия после потери устойчивости. Для случая шарнирного опирания краев (18) для уравнения (16) можно 
записать точные решения: 

(56° +27? +с.п“ ) — дп? +К 


: 1}. 1 
м о =зшлш| х-- |зшли| у-- |, 21 
‚вит = (5-1), в 


Ро," = пт? 


где т и п произвольные натуральные числа. 

Из анализа формул (21) следует, что при различных отношениях длин краев пластины 6, физических 
параметров ортотропности материала пластины с., с›, параметра 4 сжимающей нагрузки по оси у, и 
коэффициента постели К (характеризующего реакцию основания) критическому значению нагрузки 


р» = п Ри может соответствовать одна или несколько полуволн в форме равновесия (21). При этом 
туп > ь 


возможно так называемое выпучивание по двум собственным формам, когда, например, 


р. = три, =р. =р»,. Ниже приводятся такие случаи. Для решения задачи (16) с краевыми условиями 


т‚п >0 
свободного защемления краев (17) применяется вариационный метод в сочетании с конечно-разностным 
методом, который был обоснован в работах [9, 10] и легко переносится на случай ортотропной пластины. При 
этом разработанная для этой задачи компьютерная программа тестировалась на решении задачи (16), (18). 
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Кроме этого, разработана компьютерная программа решения вариационным методом в сочетании с конечно- 
разностным методом задачи (16) с каждым из краевых условий (19), (20). Для этих случаев были 
сконструированы тестовые примеры соответствующих краевых задач. 

В работе [7] в задаче об асимптотическом интегрировании уравнения малых колебаний длинной (р << а) 


прямоугольной ортотропной пластины построено разложение решения по степеням малого параметра 6 =Ъ/а и 
показано, что главный член разложения определяется из уравнения колебаний балки с модулем упругости 
Е=Е,. Повторяя рассуждения работ [7] и [10] будем искать решение задачи (16) с краевыми условиями (19) 


или (20) в виде рядов: 


со со 


р= У рыб? .м= У мб” . (22) 


1=0 1=0 
Можно показать, что главный член асимптотического разложения (22) решения задачи (16) с краевыми 
условиями (19) или (20) определяется из задачи в безразмерных переменных о критических нагрузках сжатой 
упругой балки на упругом основании со свободно защемленными или подвижно шарнирно опертыми краями: 


(1-и?) му К, + рожи =0, (23) 
\=\, =0, (24) 
У, = УМ =0. (25) 


Здесь штрихом обозначена обыкновенная производная функции, зависящей от переменной х. Можно построить 
алгоритм определения всех коэффициентов слагаемых разложений (22), но авторы не будут этого делать в 
данной работе. 

При переходе к размерным переменным из (23)-(25) получаем задачу о критических нагрузках сжатой балки 


с модулем упругости Е =Е,. Задача на собственные значения (23) с краевыми условиями (25) имеет точное 


решение: 
в. лит м, (х) =зт(лшх). (26) 


В случае краевых условий (24) задача (23) решается вариационно-разностным методом. 
Результаты исследования. Проведены численные расчеты критических значений параметра нагрузки р, 
сжимающей вдоль направления оси х при фиксированных значениях коэффициента постели К, и 
фиксированных значениях параметра нагрузки 4, сжимающей вдоль направления оси у и заданных значениях 
параметров ортотропии материала пластины. Построены соответствующие критическим нагрузкам формы 

равновесия пластины. 
Пример 1. Рассмотрим пример квадратной ортотропной пластины из работы [16] «фанера первого типа» — 
трехслойная фанера: если направить ось х вдоль волокон рубашки (верхнего и нижнего из трех слоёв), то 
1,4 


модули Юнга вдоль волокон и поперек соответственно равны Е, =1,4-10°кг/см°и Е, = 15 `10’кг/ см? ; 

коэффициенты Пуассона — у =0,46 и у, =0,46/12; модуль сдвига для главных направлений упругости 

С =12.10’кг/см?. Тогда для задачи в безразмерных переменных (16) получим 5 Ь 1 © В 1 0,083; 
а 


О, Е 


с, - 63 — 0,207 учитывая приведенные выше формулы перехода к безразмерным переменным. Пусть края 


1 


пластины шарнирно оперты. Тогда из (21) при К = 0, 4 = 8,193 следует р. = п р. =р,, =р,> =6,580. Этому 


| Ти: 1 
критическому значению отвечают две формы равновесия У 1 Е ЗШЛ [х — 5 зшлд| у- 7 и 


У ах од и 
12 5 у 5 


В таблице 1 приведены значения критической нагрузки р’ потери устойчивости ортотропной пластины с 
защемленными краями, сжатой по оси х, для различных значений коэффициента постели К и различных 
значениях сжимающей нагрузки по оси у. Таблица 1 показывает, что при повышении значения коэффициента 
постели критическая нагрузка потери устойчивости сжатой ортотропной пластины также повышается. 
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Таблица 1 
Зависимость критической нагрузки р’от параметра основания К при различных 4 
Р. 

К а=0 а = 10 а=20 

0 47,455 38,011 5,298 
1000 112,142 105,841 88,1469 
2000 146,99 144,252 122,039 
3000 164,809 163,198 141,578 
4000 177,741 176.455 175,106 
5000 189,002 187,843 186,65 


Пример 2. Рассмотрим ту же ортотропную пластину из примера | и направим ось х поперек волокон 


1,4 


рубашки. Тогда модули Юнга вдоль осей х и у соответственно равны Е, = т -10°кг/см?; В, =1,4.10°кг/см”; 


коэффициенты Пуассона — у, =0,46/12 и у, =0,46; модуль сдвига для главных направлений упругости 
С=1,2-10’кг/см’. Тогда для задачи в безразмерных переменных (12) получим И: с, = В& —12; 
а р Е 


с, и =2, 481. Пусть выполняются краевые условия (13) — края пластины свободно защемлены. В результате 
1 

численных расчетов с помощью компьютерной программы решения задачи (12), (13) для данной ортотропной 

пластины построены зависимости между нагрузкой 4 по осиу и критическим значением р, нагрузки р по оси х 


при фиксированных значениях остальных параметров этой задачи. В таблице 2 приведены значения (Ч,р +) при 


К=0. 
Таблица 2 


Зависимость между аи р* прик=0 


№ пп 1 2 3 4 5 6 
а 0 100 200 300 400 500 
В 254,388 | 241,867 | 227/789 | 208,778 | 166,462 77,105 


Отметим, что всем шести точкам (9,р*) отвечает форма равновесия, постепенно трансформирующаяся от 
вида, представленного на рис. 1 а, до вида на рис. 1 в. На рис. 1 проиллюстрированы первая, третья и шестая 
собственные функции, соответствующие критическим точкам (9,р*) таблицы 1. Эти собственные функции 


показывают формы равновесия пластины после потери устойчивости (формы потери устойчивости). 


Рис. 1. Формы потери устойчивости пластины: а — при 4 = 0, р = 254,388; 
б — приа = 200, р = 227,789; в — при 4 = 400, р = 166,462 


Пешхоев И. М. и др. Критические нагрузки равномерно сжатой ортотропной прямоугольной пластины на упругом основании 


В результате численных расчетов с помощью компьютерной программы решения краевой задачи (16), (17) 
для данной ортотропной пластины найдены значения параметров К = 200; 4=0 при которых собственному 
значению р*=255,6 отвечают две собственные функции \, и \,. Этот случай соответствует потере 


устойчивости пластины по двум собственным формам. На рис. 2 представлены графики соответствующих 
этому случаю форм потери устойчивости. 


Рис. 2. Две формы потери устойчивости: а — \,; 6 — м, 


Пример 3. Рассмотрим ортотропную пластину с параметрами Е, =1,4.10*кг/см”. Е, о 6 бела 
- 12 
коэффициенты Пуассона — у =0,46 и у, =0,46/12; модуль сдвига для главных направлений упругости 
а Ь 
С=Ь2.-10?кг/см”. Пусть длина пластины а=5, ширина 6 =1, края х= = — защемлены, а края у= 5 == 


1 
свободны, параметр 5=5 (пластину можно считать длинной), коэффициент постели К=350. Численный 


расчет решения задачи (16) с краевыми условиями (17) дает значение критической нагрузки р=37,42 ‚ график 


соответствующей собственной функции представлен на рис. 3 а. Численный расчет решения задачи (23), (24) 
при у, =0,46 и К=350 для балки длиной а=5 дает р, =39,01, при этом собственная функция имеет вид, 


представленный на рис. 3 6. Как видим, критическая нагрузка потери устойчивости пластины совпадает с 
критической нагрузкой балки с погрешностью |Р-Ро| =|37,42-39,01] =0,04. Собственные функции задач (16), 


(17) и (23), (24) имеют одинаковое число полуволн. 
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Рис. 3. Формы потери устойчивости: а — \ (х, у) ‚ сжатой ортотропной длинной пластины; б — \у (х) ‚ сжатой балки 


Обсуждение и заключения. Исследована задача о критических нагрузках сжатой в двух направлениях 
ортотропной пластины, лежащей на упругом основании. Задача сведена к решению параметрической линейной 
краевой задачи на собственные значения. В случае краевых условий подвижного шарнирного опирания 
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приведены точные формулы собственных значений и собственных функций, а в случае свободного защемления 
краев применен вариационный метод в сочетании с конечно-разностным методом и построена компьютерная 
программа решения задачи. Получены следующие выводы: 

1. При повышении составляющей сжимающей нагрузки 4 вдоль оси у снижается величина критического 
значения р * нагрузки, сжимающей пластину вдоль оси х. 


2. Если ортотропная пластина сжимается нагрузкой вдоль направления, которое соответствует большей 
изгибной жесткости, то критическое значение потери устойчивости больше, чем критическое значение 
действующей вдоль направления меньшей изгибной жесткости сжимающей нагрузки. 

3. При наличии упругого основания большему значению коэффициента постели К > 0 соответствует большее 
значение критического значения сжимающей нагрузки. Наличие упругого основания повышает несущую 
способность сжатой пластины. 

4. В случае длинной ортотропной пластины на упругом основании критическую нагрузку сжатой пластины 
можно вычислять приближенно как критическую нагрузку сжатой балки на упругом основании с модулем 
упругости, равным модулю упругости пластины в продольном направлении. 
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